Regularities and symmetries of subsets of collective 0"*" states 
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The energies of subsets of excited 0"*" states in geometric collective models are investigated and 
found to exhibit intriguing regularities. In models with an infinite square well potential, it is found 
that a single formula, dependent on only the number of dimensions, describes a subset of 0"*" states. 
The same behavior of a subset of 0^ states is seen in the large boson number limit of the Interacting 
Boson Approximation (IBA) model near the critical point of a first order phase transition, in contrast 
to the fact that these 0^ state energies exhibit a harmonic behavior in all three limiting symmetries 
of the IBA. Finally, the observed regularities in 0"*" energies are analyzed in terms of the underlying 
group theoretical framework of the difi'erent models. 



I. INTRODUCTION 

One of the overarching themes of the science of com- 
plex many-body systems is to understand the remarkable 
regularities they often exhibit and try to relate these to 
underlying symmetries of the system. In nuclei, this chal- 
lenge is approached through the use of geometric and al- 
gebraic models that describe collective behavior of the 
nuclear system. There are a large number and variety of 
such models, each with seemingly unique properties and 
predictions. Nevertheless, careful analysis often shows re- 
lations among such models that have escaped notice and 
therefore leads to a better understanding of their mutual 
interrelationships and often, to experimental tests and 
constraints on their applicability. 

In the present work, we will focus on the properties of 
subsets of 0+ states in nuclei. In general, 0+ states are 
of fundamental importance since they are easily observed 
experimentally in reactions such as few nucleon transfer 
S| or beta decay [3| . Although not all 0+ states are 
collective in nature, they are always intrinsic excitations 
of the ground state condensate, and are free of some of 
the complications (such as centrifugal effects) present in 
other states. In the present work, we focus on the prop- 
erties of a subset of collective 0+ states. We will show 
that broad classes of seemingly diverse models actually 
yield identical predictions for that subset of O"*" states. 
One upshot of this analysis will be the development of 
extremely simple, analytic eigenvalue expressions which, 
in one case, depend only on the dimensionality of the sys- 
tem and, in another, turn out to transcend the symmetry 
structure. 

Our approach primarily exploits two classes of mod- 
els, namely the Interacting Boson Approximation (IBA) 
model and geometric descriptions of nuclei at critical 
points of quantum phase transitions in their equilibrium 
structure. Therefore, we start by briefly recalling per- 
tinent features of these models. The Interacting Bo- 
son Approximation model [3] describes collective struc- 
ture in terms of bosons of angular momentum zero (s- 
bosons) and two (d-bosons) in the framework of an over- 
all U(6) symmetry. Emerging from the U(6) symmetry 
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FIG. I: (Top) Symmetry triangle with the three IBA dynam- 
ical symmetries placed at the vertices. The phase transition 
region is indicated by the slanted lines. (Bottom) Similar tri- 
angle in a geometrical framework. The critical point symme- 
tries E(5) and X(5) are located close to the phase transition 
region of the IBA. The models E(5)-/3^''' and X(5)-/3^'= span 
structures between a vibrator and the critical point symme- 
tries. 



group structure are three dynamical symmetries which 
have long been benchmark paradigms of structure : U(5), 
which gives vibrational structure characteristic of spher- 
ical nuclei, SU(3), which describes axially symmetric de- 
formed rotors, and 0(6), which pertains to deformed nu- 
clei that are soft with respect to axial asymmetry (7- 
soft). Shape/phase transitions in atomic nuclei were dis- 
cussed [1] many years ago in the classical analog 0, Q of 
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FIG. 2: Experimental [E{0+)-E{6+)]/ E{2+) values plotted as 
a function of the -R4/2 = Eji^) / E(2t) ratio for all even-even 
nuclei with Z > 40. Data from Ref. fs^]. 



the IBA. 

To visualize these limiting symmetries and the tran- 
sitions between them, it is common to place them at 
the corners of a symmetry triangle [9^], as shown in Fig. 
l(top). In the IBA framework, a point of first order phase 
transition occurs between U(5) and SU(3), while a point 
of second order phase transition occurs between U(5) and 
0(6). The triangle is divided into two regions, spherical 
and deformed, by a narrow shape coexistence region [To| 
extending around the line of first order phase transition 
connecting the two points mentioned above. In the classi- 
cal limit of the IBA, obtained through use of the intrinsic 
state formalism 0, Q , one can use (lT| Landau theory to 
delineate a similar phase transitional behavior. 

More recently, phase transitions have been investigated 
in a geometrical framework. The critical point symme- 
tries E(5) [13] and X(5) [l^l have been developed to de- 
scribe phase/shape transitions between vibrational to 7- 
soft and vibrational to axially symmetric deformed, re- 
spectively using special solutions of the Bohr Hamil- 
tonian |14l |. These solutions are analytic and parame- 
ter free, except for scale. In Fig. 1 (bottom), we indi- 
cate their position close to the critical point of a phase 
transition in a symmetry triangle for geometrical mod- 
els. The concept of critical point symmetries (CPS) is 
supported by the observation of nuclei exhibiting such 
properties [ll [11, [13, [M H H^l- Their success has 
spawned the development of numerous additional geo- 
metrical models, several of which offer analytic solutions 
and cover a wider range of structures both before and 
after the phase transitional point. 

Links between the geometrical approach and the IBA 
formalism have also found renewed interest. A power- 
ful method for solving the Bohr Hamiltonian numeri- 
cally has recently been developed [l^], leading to an al- 
gebraic collective model [l^ . Examples of the use of this 
method have recently been presented (23j . The relation- 
ship between the algebraic collective model and the dif- 



ferent limiting symmetries of the IBA has been studied 
in Refs. [23, [l^- The connection between geometrical 
models spanning structure near E(5) and the IBA has 
also been investigated [2^. 

A key prediction of the CPS involves the energy of the 
first excited 0+ state. The nature of low-lying O"*" states 
is critical to understanding the structure of nuclei and 
changes in structure 3j27l. [28j . While there is some de- 
bate as to the nature [23| of low-lying O"*" states, the fact 
that the energies of 0+ states evolve rather smoothly as 
a function of changing structure cannot be ignored. This 
has been pointed out previously in Ref. [13] ■ We illustrate 
this in a similar way in Fig. 2, plotting the relationship 
between a level of the ground state band, E{6^), and 
the first excited 0"^ state, i5(0j), as a function of i?4/2 = 
E{'i^)/E{2f) for aU even-even nuclei with Z > 40. De- 
spite the enormous range of structures encompassed in 
the plot, an overall compact trajectory emerges. 

The primary purpose of the present work is to inves- 
tigate the energies of 0"*" states in a wide range of mod- 
els. In some cases the analysis applies to all collective 
0+ states, in others to classes of 0+ states that act as 
bandheads for major families of states. For a broad class 
of models, we will discuss some remarkable regularities, 
develop analytic expressions for the eigenvalues of these 
states that depend only on the dimensionality of the sys- 
tem, relate these results to more general models, and 
discuss the group theoretical properties underlying these 
regularities. We begin within the framework of the differ- 
ent solutions of the Bohr Hamiltonian. We then do a sim- 
ilar analysis within the framework of the IBA and then 
finally investigate the links between these two different 
approaches. Some of this material has been previously 
summarized 1321 l33l . 



II. 



0+ STATES IN SOLUTIONS OF THE BOHR 
HAMILTONIAN 



Numerous models are emerging which provide a rea- 
sonable description of nuclei using an infinite square well 
potential. E(5) and X(5) are special solutions of the Bohr 
Hamiltonian describing collective nuclear properties in 
terms of the shape variables /3 and 7. Both take the po- 
tential in P as an infinite square well but use different 
potentials in 7; X(5) uses a harmonic oscillator potential 
in 7 which has a minimum at 7 = 0° whereas E(5) takes 
a potential independent of 7. Additional solutions which 
make use of infinite square well potentials in (3 include 
Z(5) [13], which uses a harmonic oscillator potential in 7 
with a minimum at 7 = 30°, Z(4) (ssj where 7 is frozen 
to 30°, and X(3) [sgI] where 7 is fixed at 0°. 

In each of the infinite square well solutions, the energy 
eigenvalues are proportional to the squares of roots of 
the Bessel functions, J,^{z), where the order v is differ- 
ent for each solution. The orders of the Bessel functions 
obtained in the E(5), X(5), Z(5), Z(4), and X(3) models 
are summarized in Table I, along with the dimension, D, 
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TABLE I: Order v, dimension, D, of the model space and 
V for = 0"^ states in the geometrical models E(5), X(5), 
Z(5), Z(4), and X(3). J is the spin of the level, r = J/2 
(in the ground state band), and n„ is the wobbling quantum 
number [33, which is zero for 0"*" states. 



Model 








!/(J=0+) 


E(5) 






5 


3 
2 


X(5) 




9 
4 


5 


3 
2 


Z(5) 


^J(J+4)+3n„(2J- 
2 


-n„) + 9 


5 


3 
2 


Z(4) 


-^J(J+4)+3n„(2J- 




4 


1 


2 




X(3) 




1 
4 


3 


1 
2 



of each model and the value of v for J'^ = 0+ states. The 
dimension effectively refers to the number of degrees of 
freedom of the model. For example, the five dimensional 
models are described by /3, 7 and the three Euler angles. 

In the X(3) and Z(4) models, we consider all excited 
0+ states. In X(5) and Z(5), the solutions are obtained 
through an approximate separation of the /3 and 7 de- 
grees of freedom. We consider those 0+ states arising 
from the (3 solution, as they are directly related to the 
infinite square well potential. In the E(5) solution, we 
consider those 0+ states with t = 0, that is, those 0+ 
states which correspond to base states on which major 
families of levels are built. 

Traditionally, the excitation energy of the first excited 
2+ state is used to set the overall scale for these models. 
However, in some cases, using a different normalization 
can reveal physics not otherwise very evident. In par- 
ticular, one can sometimes see relations among states of 
the same angular momentum by normalizing to the first 
excited state of that spin. Hence, here, we scale to the 
energy of the first excited 0+ state, 0^ . It turns out that 
this approach allows the relative energies of these excited 
0+ states to be well described by simple formulas. For 
the 0+ states we use the usual notation, 0^, where m ~ 1 
corresponds to the ground state, m = 2 denotes the first 
excited 0+ state, and so on, within each of the subsets of 
0+ states, described above. 

The energies of the 0+ states in the X(3) model are 
given in Table II. Normalizing to the first excited 0+ 
state, O2 , the energies are described exactly by 



E{0+)=An{n + 2), 



1, 



(1) 



where n gives the sequencing of 0+ states defined such 
that the first excited 0+ state corresponds to n = 1, and 
A is a scaling factor. 

In the Z(4) model, normalizing to the first excited 0+ 
state, the 0^ states increase approximately as 



TABLE II: Energies of 0^ states in different geometrical mod- 
els using an infinite square well potential. Columns labelled 
2^ (Oj) are normalized to the first excited 2"'" (0^) energies. 
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ot 


2.87 


1.0 


2.95 


1.0 


3.03 


3.91 


5.65 


1.0 


ot 


7.65 


2.67 


7.60 


2.57 


7.58 


9.78 


14.12 


2.50 


ot 


14.34 


5.00 


13.93 


4.71 


13.64 


17.61 


25.41 


4.50 


ot 


22.95 


8.00 


21.95 


7.43 


21.22 


27.39 


39.53 


7.00 


ot 


33.47 


11.67 


31.65 


10.72 


30.31 


39.12 


56.47 


10.00 



where again, A is a scaling factor dependent on the par- 
ticular model. 

The energies of 0+ states in the E(5), Z(5), and X(5) 
models, normalized to the 2^ state, are given in Table II. 
While these energies are quite different, by normalizing 
to the first excited 0"^ state, the models produce exactly 
identical results, as seen on the right of Table II. These 
energies very closely follow the simple formula 



E{0+) = An{n + 3), 



n = m — 1. 



(3) 



Equation ^ is not an exact description of the 0+ en- 
ergies in E(5), Z(5), and X(5); however, it does provide a 
very accurate approximation. Through n=10, the model 
0"*" energies deviate from the expression given in Eq. ([3]) 
on the order of less than 0.1 %. 

The above empirical results can all be combined into 
a single, simple formula describing the 0^ states in any 
model with an infinite square well potential by 



i?(0+) = AnU + 



D + 1 



1 



(4) 



where D is the number of dimensions and A again de- 
pends on the model. Eq. ([4]) is exact only for D — 3. 
As mentioned previously, for other values of D and low 
values of i^, it is a very accurate approximation. 

Equation (4) stems from particular relations between 
the zeros of the Bessel functions which are involved in 
the solutions to the infinite square well models. Given 
the Bessel function J„{z), with roots Zj, s = 1, 2, 3, 
... we empiricaly observe that the following approximate 
relation holds 



(5) 



E{0+) ^ An{n + 2.5), n = m-l, 



(2) 



where n = s — 1. We call E^{n) the spectrum of the roots 
of J^, since this quantity corresponds to energy spectra 
in models describing atomic nuclei. 

The relation given in Eq. (5) is exact only in the case 
1/ = 1/2, as one can see [38| from the expansions of roots 
of Bessel functions given in [s^ (Eq. 9.5.12). Numer- 
ical results are shown in Table HI. It is clear that the 
approximation deteriorates rather slowly with increasing 
n (while keeping v constant), while it deteriorates faster 
with increasing v (while keeping n constant). 
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TABLE III: Exact spectra of several Bessel functions Jv (labeled by E,j{n)) compared to the corresponding n(n+3;) approximate 
expressions (labeled by x). J3/2 occurs in E(5), X(5), Z(5). J1/2 occurs in X(3). Ji occurs in Z(4). Jo, Ji, J2, J3 occur in the 
pairing case [40| |. 



n 


Eo{n) 


x=1.5 


1 


1.000 


1.0 


2 


2.799 


2.8 


3 


5.398 


5.4 


4 


8.796 


8.8 


5 


12.993 


13.0 


6 


17.990 


18.0 


7 


23.787 


23.8 


8 


30.383 


30.4 


9 


37.779 


37.8 


10 


45.974 


46.0 



gi/2(n) x^2 



I. 000 
2.667 
5.000 
8.000 

II. 667 
16.000 
21.000 
26.667 
33.000 
40.000 



I. 000 
2.667 
5.000 
8.000 

II. 667 
16.000 
21.000 
26.667 
33.000 
40.000 



=2.5 



1.000 
2.572 
4.715 
7.430 
10.716 
14.574 
19.004 
24.005 
29.577 
35.721 



1.000 

2.571 

4.714 

7.429 

10.714 

14.571 

19.000 

24.000 

29.571 

35.714 



g3/2(n) x^3 



1.000 
2.500 
4.499 
6.999 
9.998 
13.497 
17.496 
21.995 
26.993 
32.492 



1.0 
2.5 
4.5 
7.0 
10.0 
13.5 
17.5 
22.0 
27.0 
32.5 



E.2{n) 



=3.5 



1.000 
2.443 
4.329 
6.659 
9.433 
12.651 
16.313 
20.418 
24.967 
29.960 



1.000 
2.444 
4.333 
6.667 
9.444 
12.667 
16.333 
20.444 
25.000 
30.000 



Er,/2{n) x=4 



I. 000 
2.397 
4.192 
6.385 
8.977 

II. 968 
15.357 
19.145 
23.332 
27.918 



1.0 
2.4 
4.2 
6.4 
9.0 
12.0 
15.4 
19.2 
23.4 
28.0 



Esin) x=4.5 



I. 000 
2.358 
4.077 
6.157 
8.599 

II. 403 
14.568 
18.095 
21.983 
26.233 



I. 000 
2.364 
4.091 
6.182 
8.636 

II. 455 
14.636 
18.182 
22.091 
26.364 



For the Bessel function Y^{z), with roots Zg, s = 1, 
2, 3, ... (not used in this paper), we observe that the 
following approximate relation holds 

^•^W - 72 3 —3 ' (6) 

where n = s — 1. This formula is exact only in the 
case ly — 1/2. Similar to the results obtained for the 
Bessel function J,y{z) relations, the approximation dete- 
riorates rather slowly with increasing n (while keeping v 
constant), while it deteriorates faster with increasing v 
(while keeping n constant). 

These results are applicable outside of the models dis- 
cussed above. As an example, Eq. ([4]) applies to a recent 
model j40| describing the critical point of a pairing vibra- 
tion to pairing rotation phase transition. Here the ener- 
gies of 0+ states, which span two degrees of freedom, the 
excitation energies of a particular nucleus and the masses 
along a series of even-even nuclei, can be described. In 
addition, these results would be applicable to hadronic 
spectra, which have recently been described [iH in terms 
of roots of Bessel functions. 

We can go further and generalize Eq. ^ for a broader 
range of models. The Bohr Hamiltonian can also be 
solved with potentials in (3 of the form V ~ /S^'^, giv- 
ing the so-called E(5)-/32'= model [H and the X(5)-/32'= 
model jisl, using the 7 dependence characteristic of E(5) 
and X(5), respectively. These models allow for a descrip- 
tion of structure between vibrational- like and the infinite 
square well solutions by increasing the power of /? in the 
potential. For example, in the E(5)-/3^'' case, gives 
the vibrational limit and as the power of f3 goes to infin- 
ity, the E(5) solution is reached. The evolution of both 
models is included schematically in Fig. 1 (bottom). The 
predicted 0+ energies of the (3'^'^ models are plotted in 
Fig. 3(top) normalized to the first 2+ state energy. As 
evident from Fig. 3 (top), with increasing powers of (3 in 
the potential, the 0"*" energies evolve gradually towards 
the infinite square well predictions. However, again, the 
E(5) and X(5) related models seemingly give different 
results. If instead, we normalize each energy to that of 
the first excited 0"*" energy, these models produce exactly 



identical results for a given /3 potential, as seen in Fig. 
3 (bottom). The normalized E{0^) energies can be re- 
produced with a generalized version of Eq. ^ given by 

E{0+) = An{n + x), n = m - 1, (7) 

where x is some number. The values of x obtained by 
fitting the first two 0+ state energies in each model with 
Eq. ([7|) are included in Fig. 3. For a harmonic oscillator 
potential in (3, x — 00, since the bandhead 0^ energies 
increase linearly (i.e., when considering the ratio of ener- 
gies the term in parenthesis in Eq. (7) disappears). As 
the power of P in the potential is increased, the value of x 
decreases, reaching the limiting value of 3 for the infinite 
square well. 

III. 0+ STATES IN THE EUCLIDEAN 
ALGEBRAS E(n) 

In the solutions of the Bohr Hamiltonian with an infi- 
nite square well potential in the (3 degree of freedom, the 
regularities observed for 0+ states can be related to the 
second order Casimir operator of E(D), the Euclidean 
group in D dimensions. In order to see this, one can con- 
sider in general the Euclidean algebra in D dimensions, 
E(Z3), which is the semidirect sum [l^l of the algebra 
T/) of translations in D dimensions, generated by the 
momenta 



and the SO(D) algebra of rotations in D dimensions, gen- 
erated by the angular momenta 

f d d \ 

symbolically written as E(D) = Td (Bs SO(D) [ist . 

The generators of E(D) satisfy the commutation rela- 
tions 

[P„ P,] = 0, [P„ Ljk] = i{S^kPJ - kjPk). (10) 
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FIG. 3: (Color online) (Top) Excited 0^ state energies in the 
E(5)-/3^'' and X(5)-/3^* models, normalized to the 2'1 state en- 
ergy. (Bottom) Same as top, but with the energies normalized 
to the first excited 0"*" state energy. IW is the infinite square 
well potential used in the original E(5) and X(5) solutions. 



[Lij,Lki] = i{SikLji + SjiLik — SiiLjk — SjkLu). (11) 



the eigenf unctions of which are the Bessel functions 
/(r) = Mkr) m. 

The "radial" equations in the infinit e sq uare well mod- 
els E(5) [Tl, X(5) Z(5) [Sl, Z(4) and X(3) [s^ 
are obtained, after the transformation of Eq. (fTH)) has 
been performed, in the form of Eq. p^ . with the order 
v summarized in Table I. 

In E(5), Eq. (fT4|) and the corresponding order ly given 
in Table I coincide with D = 5 and lu — t, where t{t + 3) 
represents the eigenvalues of the second order Casimir 
operator of SO (5). Thus all states obey Eq. p^ . 

In X(5), where again D — 5, Eq. (fT4|) and the cor- 
responding order given in Table I would agree for 
J(J + l)/3 — Lj{uj + 3). This does not hold for any J 
in general, but it is satisfied for J = = oj. Thus all 0+ 
bandheads obey Eq. (|15p . 

In Z(5), where again D — 5, Eq. (|14p and the corre- 
sponding order ly given in Table I would agree for = 
and for J(J-|-4)/4 = uj{u! + 3). Again, this does not hold 
for any J in general, but it is satisfied for J = = cj. 
Thus 0+ bandheads with n^, = obey Eq. P3|) . 

In the case of Z(4) [Hi, in which L» = 4, Eq. ^ and 
the corresponding order v given in Table I for Uw = 
obtain the form i/ = oj + 1 and v = J/2 + 1, respectively. 
They agree for J — 2lu, as already known [S^j; therefore 
states with any even J and rt^ = obey Eq. (|15p . 

In X(3), where D — 3, Eq. fH]) and the corresponding 
order v given in Table I would agree for J(J + l)/3 = 
uj{ll! + 1). Once again, this does not hold for any J in 
general, but it is satisfied for J — — uj. Thus 0^ 
bandheads obey Eq. (fT5|l. 

The above situation is similar to a partial dynamical 
symmetry [47] of Type I [48] , where some of the states 
(the 0"*" states in the present case) preserve all the rele- 
vant symmetry. 



From these commutation relations, the square of the total 
momentum, P^, is a second order Casimir operator of the 
algebra, while the eigenfunctions of this operator satisfy 
the equation 



1 d D-id lu{lo + D-2) 



,,D-1 



dr 



dr 



F{r) = k^F{r), 
(12) 

where on the left hand side of Eq. the eigenvalues 
of the Casimir operator of SO(D), uj{uj + D — 2) appear 
[46|]. Using the transformation 



F(r) 



,,(2--D)/2 



fir), 



(13) 



IV. THE IBA HAMILTONIAN AND 
SYMMETRY TRIANGLE 

In order to describe a wider range of structures, it is 
useful to use a more general collective model than the 
specific solutions described above. To this end, we exploit 
the IBA model, which covers a gamut of structures with 
an economy of parameters. To do so, we use an IBA 
Hamiltonian of the form [491 



(1 - C)nd - 



C 



4iVB 



(16) 



and 



^ UJ + 



D 



Eq. (fT2|) can be written as 



— + -— + e--]f{r)^0, 



(14) 



(15) 



where fid ^ ■ d, = (s^d + d^s) + x(d^d)(2), 
the number of valence bosons, and c is a scaling factor. 
The above Hamiltonian contains two parameters, ^ and 
X, with the parameter ^ ranging from to 1, and the 
parameter x ranging from to — \/7/2. The U(5) sym- 
metry is given by C = 0, any x, the SU(3) symmetry 
hy ( = 1 and x = — v^/S, and the 0(6) symmetry by 
C = 1 and X = 0. With this parameterization, the entire 
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symmetry triangle, shown in Fig. 1, can be described, 
along with each of the three dynamical symmetry limits. 
Calculations in this work have been performed with the 
code IB AR [s^, IMl , which has recently been developed 
to handle large boson numbers. 

In Section V, we discuss the fact that, in all three lim- 
iting symmetries of the IBA, the energies of certain sub- 
sets of 0+ states exhibit harmonic behavior in the limit 
of large boson numbers. (This result has also been de- 
rived [S^lsst using the coherent state formalism). In con- 
trast, we will show, in Section VI, that near the critical 
point, the same subsets of 0+ states exhibit the n{n + 3) 
behavior found in the framework of geometrical models 
in Section II. Furthermore, near the critical point these 
0+ states in the large boson number limit of IBA ex- 
hibit certain degeneracies with alternate members of the 
ground state band, calling for further investigation. 



V. 0+ STATES IN THE LIMITING 
SYMMETRIES OF THE INTERACTING BOSON 
MODEL 

We begin an analysis in the IBA framework ^5] by look- 
ing at the three dynamical symmetry limits, and analyz- 
ing the behavior of the 0+ states in the analytic formu- 
lae appropriate to each, especially in the large Nb limit. 
Again, we consider a particular subset of 0+ states, look- 
ing for simple patterns common to all three symmetries, 
despite the diversity of the structures they describe. 

In the case of U(5), states are labeled by their quan- 
tum numbers v and riA, where v is the seniority and rtA 
is the number of triplets of bosons coupled to angular mo- 
mentum zero. There are two classes of O'^ states, those 
with riA = 0, and those with riA — 1,2,3. . . which are 
always found degenerate with 3^ states. In the present 
work, we consider those states with riA = 0, that is those 
0+ states that are not degenerate with 3+ states. These 
states correspond to base states on which major families 
of levels are built. In the U(5) limit, the energies of the 
0+ states with riA=0 arc proportional to the number of 
d bosons, n^, corresponding to their respective phonon 
number (terms proportional to are also allowed, but 
are omitted in the present consideration); thus, the en- 
ergies increase linearly. 

In the SU(3) hmit of the IBA, the position of the 0+ 
bandheads is determined by the second order Casimir op- 
erator of SU(3). The eigenvalue expression for 0+ states, 
in terms of the representation labels (A, /i), is given by E 
= a[A^ + fM^ + Xfi + 3(A -I- /i)]. Here, we consider all 0+ 
states. Taking the O"*" state which belongs to the {2Nb, 0) 
irreducible representation (irrep) at zero energy, and nor- 
malizing to the first excited 0+ state, which belongs to 
the {2Nb — 4, 2) irrep, we find for the lowest O"*" states 
the results shown in Table IV. From Table IV it is clear 
that at large boson numbers Nb, we have two states with 
normalized energy 2, two states with normalized energy 
3, three states with normalized energy 4, and so on. In 



other words, for large Nb, the energies of the O"*" states 
in the SU(3) limit of the IBA grow linearly. 

In the case of 0(6), states are labelled by their quan- 
tum numbers a and r. One set of excited 0''' states is 
found within the multiplet structure of a given a family, 
has r values of 3 or larger and always appears degener- 
ate with J = 6"^, 4+, and 3+ states. The other set of 
excited 0+ states forms the bandheads of the different a 
families, has r = and are not degenerate with other 
states in the spectrum. We consider only those states 
with T—Q. In the 0(6) limit of the IBA, the positions 
of the 0"'' bandheads are determined by the second order 
Casimir operator of SO (6). The eigenvalue expression for 
T ~ 0, J ~ states in terms of the major family quan- 
tum number a is E = a a{a + 4), with a = Nb, Nb — 2, 
. . . , or 1. Taking the 0+ state which belongs to the 
{Nb) irrep at zero energy, and normalizing to the first 
excited 0+ state belonging to the {Nb — 2) irrep, we ob- 
tain the results shown in Table IV. We observe that O'^ 
bandheads in the 0(6) limit of the IBA also grow linearly 
in the limit of large boson numbers Nb ■ 

For all three IBA dynamical symmetry limits, the en- 
ergies of these sets of 0+ states are given by 

E = An (17) 

in the large Nb limit. Thus, a single simple formula 
applies to all three dynamical symmetry limits of the 
IBA despite the fact that each describes very different 
structures. It is interesting that this harmonic behavior, 
appearing as a general feature of IBA spectra, not only 
at the three limiting symmetries but also in intermediate 
situations [H, [5^ , is strongly violated near the critical 
point, as we shall see in the next section. 

VI. 0+ STATE ENERGIES AND 
DEGENERACIES IN THE SHAPE 
COEXISTENCE REGION OF THE IBA 

It has been recently observed |33] that the line de- 
scribing the degeneracy E{<o^) = i?(Oj) (where 0^ is the 
first excited 0"'' state) in the symmetry triangle of the 
IBA for large Nb {Nb = 250) falls within the coexis- 
tence region of spherical and deformed shapes, slightly 
to the right of the critical line representing the first or- 
der phase transition between U(5) and SU(3). Similar 
results are obtained for the E{lQf) = E{0+) and £^(14+) 
= E{Q^) degeneracies. These degeneracies are interest- 
ing not only because they can possibly be associated with 
underlying symmetries but also because the degeneracy 
between -E(Oj) and E{(!>1) found near the critical point 
of the IBA is also approximately given by the X(5) crit- 
ical point model. In what follows, we investigate further 
the degree to which the IBA predictions near the critical 
point are related to simple analytic formulas. 

Along the U(5)-SU(3) leg (x = -\/7/2) a de generacy 
between the 6^ state and the state occurs for C, — 
0.473 for Nb = 250. This occurs very close to, but just 
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TABLE IV: Irreducible representations (irreps) of SU(3) and 0(6) and the corresponding 0"*" bandhead energies. In the case of 
SU(3), energies are normalized to the 0"*" bandhead with (A, /i) = (2A''-4,2), while in 0(6) to the 0''' bandhead with a — {N-2). 
N stands for the boson number, Nb- 



Irrep (A,^) 


E{0+) 


Irrep (A,/^) 


SU(3) 

E{0+) 


Irrep (A,/i) 


E{0+) 


0(6) 
Irrcp(cr) E{0+) 


(27V,0) 
(2iV-4,2) 
(2Ar-8,4) 
(2Ar-12,6) 
(27V-16,8) 



1 

(4Ar-6)/(2iV-l) 
(6A^-15)/(2iV-l) 
(8Af-28)/(2iV-l) 


(2Af-6,0) 
(2iV-10,2) 
(2iV-14,4) 


(4iV-3)/(2iV-l) 
(6iV-10)/(2Af-l) 
(8iV-21)/(2Af-l) 


(2iV-12,0) 


(8iV-18)/(2Af-l) 


(TV) 
(iV-2) 1 
(iV-4) 2-{2/N) 
(iV-6) 3-{3/N) 
(iV_8) i-{8/N) 



TABLE V: Predictions of the IBA (with Nb = 250, x = 
— \/7/2, C ~ 0.473) compared to analytic expressions (see 
text). On the left, excited 0"'' energies are compared while 
on the right, energies in the ground state band are compared. 
Results are normalized to E{2f) — 1.0, the numerical factors 
accompanying n(n + 3) and J( J + 2) in the column headings 
reflecting this normalization. 





Analytic 


IBA 




Analytic 


IBA 


n 


3n(n+3) 
2 




J 


./{J+2) 
8 


EiJ) 








2 


1.00 


1.00 








4 


3.00 


3.05 


1 


6.00 


6.08 


6 


6.00 


6.08 








8 


10.00 


10.00 


2 


15.00 


14.85 


10 


15.00 


14.73 








12 


21.00 


20.23 


3 


27.00 


27.57 


14 


28.00 


26.43 








16 


36.00 


33.30 


4 


42.00 


42.55 


18 


45.00 


40.81 



beyond, the critical point {(crit = 0.472 for Nb='250) 
of the phase transition. Numerical results of the above 
IBA calculation are given in Table V and compared with 
simple analytic formulas. In the first three columns, the 
first four excited O^j states (normalized to the energy of 
the 2^ state) obtained in this calculation are compared 
to the predictions of Eq. ([3]), i.e., to the n{n + 3) formula. 
Very good agreement is obtained up to n = 4. This result 
is what might be expected given the similarity of the IBA 
coherent state energy functional at the critical point with 
an infinite square well potential. 

As mentioned above, successive ground band members 
with J > 2 and J/2 odd are nearly degenerate with 
higher lying 0"*" states. Trying to satisfy simultaneously 
the degeneracies E{6t) = E{0+) and £;(10^') = E{0+) 
with the 0"^ states obeying Eq. ([3]), and the levels of the 
ground state band obeying a general equation of the form 
E{J) ~ AJ{J + y), one obtains y = 2, i.e., the levels of 
the ground state band should grow as 

E{J) = AJ{J + 2), (18) 

where, again, A is some number. In a very crude inter- 
pretation, the J( J+2) empirical result can be thought of 



as the average of the vibrational limit, where the energies 
go as J, and the rotational limit, where the energies go 
as J(J+1). 

The relevant connection to the 0^ expression is given 

by 

J(J + 2) = 12n(n + 3). (19) 

The predictions of the above mentioned IBA calcula- 
tion are compared to the J{J + 2) predictions, which are 
normalized to the energy of the 2f state, in the right 
section of Table V. Good agreement is obtained at lower 
levels, the deviation reaching 10% at J = 18. Also visible 
in the table are the approximate degeneracies E{6i) = 
E{0+), £(10+) = EiOt), E{U+) = E{0+), E{18+) = 
E{0'^). These degeneracies hold to the 10 percent level 
for J = 18. 

In summary, IBA 0+ states (in the large boson num- 
ber limit) near the critical point on the U(5)-SU(3) line 
exhibit the same n(n + 3) behavior seen in geometrical 
models involving infinite square well potentials. Further- 
more, these 0+ states demonstrate approximate degen- 
eracies with alternate members of the ground state band, 
calling for further investigations. 

VII. CONCLUSIONS 

Working within the framework of both algebraic and 
collective models, we have investigated the energies of 
subsets of excited 0+ states, pointing out regularities 
within and similarities between the two different ap- 
proaches. For models employing an infinite square well 
potential in the f3 degree of freedom, a single formula 
is derived for a subset of excited 0"*" state energies, de- 
pendent only on the number of dimensions of the model. 
The same regular behavior for 0+ states has been found 
in IBA calculations (in the large boson number limit) 
near the critical point of the first order phase transi- 
tion between U(5) and SU(3), despite the fact that in 
all three limiting symmetries of the IBA (in the large 
boson number limit) the same 0+ states exhibit a har- 
monic behavior. Furthermore, these successive 0-|- states 
near the critical point exhibit degeneracies with alternate 
yrast states, analogous to the near-degeneracy that oc- 
curs between the first 6+ state and the first excited 0+ 
state in X(5), calling for further investigations. Finally, 
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the observed regularities in 0"*" energies are discussed in 
terms of the underlying group theoretical framework of 
the different models. 
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